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ABSTRACT 


Xne  gasaynamic  problem  or  collapsing  shocks  and  detonation  waves 
having  spherical  or  cylindrical  symmetry  is  considered  near  the  point 
or  axis  of  symmetry.  The  solution  basic  to  this  work  is  the  self -similar 
flow  of  a  collapsing  symmetrical  shock  wave  with  counterpressure  neglected. 
The  focussing  effect  as  the  flow  progresses  causes  the  front  to  accelerate 
and  its  velocity  is  singular  at  the  instant  of  collapse.  In  the  present 
work  the  perturbations ,  due  to  counterpressure  and  also  to  a  uniform  heat 
release,  which  give  rise  to  essentially  identical  mathematical  solutions, 
are  evaluated.  The  basic  self -similar  solution  is  investigated  in  detail 
over  a  range  of  values  of  the  specific  heat  ratio. 
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The  similarity  solution  to  the  problem  of  a  contracting  (implod¬ 
ing)  spherical  or  cylindrical  shock  front  propagating  into  a  uniform  gas 
at  rest  is  well  known.  As  the  shock  progresses  its  surface  area  diminishes 
causing  its  velocity  to  increase  towards  the  center  of  symmetry,  where  it 
is  infinite.  The  similarity  solution  is  valid  near  the  center  of  symmetry, 
where  the  shock  is  strong. 

\ 

In  the  present  paper  the  shock  is  replaced  by  a  contracting 
detonation  front  propagating  into  a  uniform  gas  and  releasing  a  constant 
amount  of  energy  per  unit  mass  of  gas.  At  large  distances  from  the 
center,  where  the  curvature  is  negligible,  the  detonation  is  a  Chapman - 
Jouguet  front,  l.e.,  it  travels  with  sonic  speed  relative  to  the  burnt 
gas.  The  front  accelerates  towards  the  center  of  symmetry  and  becomes 
overdriven,  the  motion  now  being  governed  more  by  the  compression  effects, 
due  to  focussing  of  the  front,  than  by  the  heat  release.  The  solution 
for  the  final  stages  is  obtained  as  a  perturbation,  of  order  the  inverse 
square  of  the  speed  of  the  front,  on  the  corresponding  similarity  solu¬ 
tion  involving  a  shock  wave.  In  the  latter  solution  the  strong  shock 
relations  are  applied  at  the  front  so  that  only  the  undisturbed  density 
enters  into  the  problem,  which  has  no  time  scale.  In  the  present  case 
of  a  detonation  the  heat  releaae  is  taken  into  account,  to  first  order 
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so  that  the  basic  similarity  hypothesis  is  unaltered,  in  the  conservation 
equations  at  the  front  Tfo  form  of  the  perturbation  sc  obtained  la 
identical  to  that  due  to  taking  into  account  the  pressure  Cor  sound 
speed)  of  the  undisturbed  gas  to  first  order.  The  disturbance  of  the 
speed  of  the  front  due  to  heat  r  lease  and  initial  pressure  are  evaluated 
for  several  spherical  and  cylindrical  cases  by  linearizing  the  equations 
of  motion.  The  solution  has  to  satisfy  the  conservation  equations  at 
the  front  and  also  be  regular  on  a  certain  characteristic.  The  basic 
and  perturbation  equations  are  integrated  numerically  by  making  use  of 
the  power  series  expansions  about  this  characteristic.  A  comparison  is 
made  with  the  results  obtained  by  the  approximate  method  as  given  by 
Whitham. 

The  results  obtained  by  Butler  for  the  Guderley  solution  are 
recomputed  and  extended.  It  is  found  necessary  to  investigate  the 
existence  and  uniquenes  of  this  solution. 

1.  INTRODUCTION 

The  unsteady  motion  of  a  perfect,  inviscid,  non-heat -conducting 
gas  is  in  general  governed  by  partial  differential  equations.  However, 
in  the  case  of  a  flow  which  is  one -dimensional,  or  spherically  or 
cylindrically  symmetric,  so  that  the  flow  variables  depend  on  a  dis¬ 
tance  coordinate  R  and  the  time  coordinate  t,  there  is  a  class  of  solu¬ 
tions  in  which  all  variables  are  functions  of  a  single  combination  of  R 
and  t,  R/ta  where  a  is  a  constant.  Such  flows  are  self -similar  (Sedov, 
1959)  and  are  governed  by  ordinary  differential  equations.  The  special 
case  a  =  i  corresponds  to  a  uniformly  expanding  or  contracting  flow,  so 
that  if  such  a  flow  is  adiabatic  then  it  is  also  homentropic,  apart 
from  entropy  jumps  across  discontinuities,  as  any  shock  wave  in  the 
flow  is  of  uniform  strength.  An  example  of  a  flow  of  this  type  is  that 
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of  a  strong  point-explosion  (Sedov,  1959;  Taylor,  1950i>)  which  involves 
aCi  expanding,  decaying  spherical  snocK  wave. 

The  problem  to  be  investigated  here  is  that  of  a  contracting 
spherical  or  cylindrical  detonation  wave  propagating  into  a  uniform 
combustible  gas.  It  is  already  known  that  there  is  no  solution  involving 
a  uniformly  contracting  front  (Selberg,  1959;  Stanyukovich,  19f0).  This 
result  will  be  deduced  later  from  investigation  of  the  integral  curves 
of  R/t,  homentropic  solutions. 

In  order  to  solve  the  problem  of  a  contracting  detonation  front 

it  will  be  necessary  to  study  Guderley' s  solution  (Guderley,  1942;  Butler,  1954) 

for  a  converging  shock  wave,  in  which  the  shock  front  accelerates  towards 

R  ■  0,  where  its  velocity  is  infinite.  If  the  shock  path  is  R  “  k  (t), 

* 

then  the  shock  speed  U  is  given  by 
*  1-1 /a 

U  -  -  \  where  0  <  a  <  1. 


The  Guderley  similarity  solution  is  valid  for  small  values  of  X,  for 
which  the  shock  is  strong  so  that  the  undisturbed  gas  pressure  can  be 
neglected.  The  flow  variables  behind  the  front  thus  depend  upon  the 
shock  speed  and  undisturbed  density  only,  which  leads  to  the  similarity 
hypothesis.  If  we  now  consider  the  effect  of  a  uniform  heat  release  in 
the  medium  as  the  front  passes  through  it,  then  this  results  in  the 
addition  of  a  finite  amount  of  energy  per  unit  mass  to  the  system  and  is 
thus  a  perturbation  on  the  Guderley  solution.  The  form  of  the  perturba¬ 
tion  can  be  deduced  as  follows.  The  particle  velocity  behind  the  shock 
is  given  by 
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Let  the  particle  velocity  behind  the  detonation  be 

*  * 

u,,  *  u  +  V 
D  s 

* 

where  V  is  supposed  small  relative  to  ug.  The  extra  kinetic  energy  per 
unit  mass,  which  is  directly  due  to  the  heat  release  and  so  must  be 
finite,  is  of  order  Vu  and  hence  the  perturbation  velocity  V  is  of  order 

3 

^  2  2+2  f  CL 

U  or  relative  to  the  basic,  shock  wave  solution.  Similarly 

the  sound  speed  perturbation  is  of  order  X  242/at  Throughout  the  flow 
in  general  the  perturbations  are  of  order  R  2+2/a.  The  effect  of  allow¬ 
ing  for  the  initial  pressure  (or  internal  energy)  of  the  undisturbed  gas 
gives  rise  to  perturbations  of  precisely  the  same  form.  Let  the  speed 
of  the  front  be  given  by 


y}-Uct  (1+0X“2+2/a) 


where  0  is  a  constant  due  either  to  heat  release  or  initial  pressure,  or 
both.  In  a  given  case  we  require  the  values  of  a,  0  to  determine  the 
path  of  the  front. 

The  evaluation  of  the  constant  parameter  a  is  performed  by 
integrating  the  equations  of  motion,  which  can  be  reduced  to  a  single 
first  order,  non-linear  differential  equation,  subject  to  certain  boundary 
conditions.  In  the  case  of  a  point -explosion  a  is  determined  simply  by 
consideration  of  the  dimensions  of  the  basic  parameters  (the  density  and 
the  energy  of  the  explosion).  However,  in  the  contracting  case  there  is 
only  one  basic  parameter,  the  density,  and  a  unique  mathematical  solu¬ 
tion  is  obtained  by  assuming  that  the  flow  is  regular  on  a  certain 
characteristic  following  behind  the  shock.  The  conservation  equations 
across  the  front  and  the  regularity  condition  on  the 
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characteristic  provide  the  two  necessary  boundary  conditions  for  the 
solution  of  the  differential  equation.  The  values  of  a  for  the  six 
cases  v  *  1.2,  1.4,  5/3,  spherical  and  cylindrical,  have  been  computed 
by  Butler.  For  an  arbitrary  choice  of  a  there  are  four  possible  solu¬ 
tions  which  satisfy  the  regularity  condition  but  in  each  case  only  one 
of  these  can  be  made  to  satisfy  the  conservation  equations  at  the  front. 
In  order  to  solve  the  equations  governing  the  perturbations  it  is 
necessary  to  recompute  the  basic  similarity  solution.  Butler's  results 
are  extended  to  y  =  3,  for  the  products  of  a  detonation,  and  in  this 
case  it  is  found  that  the  choice  of  solution  from  the  four  possibilities 
differs  from  that  for  the  lower  values  of  y.  The  changeover  from  one 
solution  to  the  other  is  discussed  in  terms  of  the  integral  curves  of 
the  system. 

To  find  the  correct  value  of  a  we  must  use  a  method  of  trial 
and  error.  The  equation  is  integrated  with  an  arbitrary  value  of  a, 
starting  at  the  characteristic  with  the  regularity  condition  satisfied, 
and  the  discrepancy  between  this  solution  and  the  shock  point  noted. 

We  repeat  the  process  with  various  values  of  a,  until  one  is  found 
passing  through  the  shock  point  to  the  required  accuracy. 

The  equations  governing  the  perturbations  are  three  simultaneous, 
linear,  first  order  differential  equations,  the  coefficients  containing 
the  basic  similarity  solution.  Again  the  solution  has  to  satisfy 
boundary  conditions  at  the  front  and  the  characteristic,  the  displace¬ 
ments  of  each  from  the  basic  solution  being  accounted  for.  However,  the 
linearity  of  the  equations  means  that  the  appropriate  solution  can  be 
evaluated  by  taking  a  certain  combination  of  any  two  linearly  independent 


solutions . 
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The  method  of  integration  of  the  basic  and  perturbation  equations 
makes  Uac  ol  the  power  series  expansions  about  the  critical  characteristic. 
In  this  way  we  avoid  any  difficulty  due  to  derivatives  being  indetermi¬ 
nate  on  direct  substitution  into  the  differential  equations.  The  solu¬ 
tion  is  developed  by  an  iterative  procedure  which  produces  approximations 
to  the  solution  and  its  derivative  in  tabular  form.  Each  iteration 
effectively  takes  into  account  another  term  in  the  power  series.  0  is 
calculated  for  y  =  1.2,  1.4,  5/3,  3  for  both  cylindrical  and  spherical 
symmetry,  and  for  heat  release  and  undisturbed  pressure.  Comparison  is 
made  with  results  obtained  by  the  approximate  method  in  the  form  given 
by  Whitham  (1958).  It  Is  known  that  this  approximate  method,  as  applied 
by  Chisnell  (1957)  in  his  "shock-area"  rule,  gives  extremely  accurate 
results  for  the  values  of  a  but  it  is  found  here  that  the  approximate 
values  of  0  by  this  method  are  much  less  accurate. 

The  equations  governing  the  motion  are  integrated  between  the 
front  and  the  characteristic,  which  is  necessary  for  the  evaluation  of 
a  and  0.  To  obtain  the  distribution  of  the  physical  variables  behind 
the  front  the  integration  would  have  to  be  continued  as  far  as  t  ■  0, 
at  which  instant  the  shock  is  at  R  ■  0  and  is  reflected.  If  all  the 
heat  energy  available  is  released  during  the  contracting  phase  of  the 
motion  then  the  front  is  reflected  as  a  shock  wave. 

Contracting  shock  waves  have  previously  been  investigated  both 

experimentally  (Perry  and  Kantrowitz,  1951)  and  numerically  (Payne, 

1957).  A  problem  having  great  similarity  to  that  of  converging  shocks 

is  that  of  cavitation  in  water,  which  has  been  studied  by  Hunter  (1960, 

1963)  and  differs  from  the  former  in  the  boundary  conditions  at  the  front 

and  the  fact  that  the  motion  is  taken  to  be  homentropic.  A  regularity 

condition  on  a  certain  characteristic  is  also  employed  to  obtain  a 
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unique  solution.  The  similarity  hypothesis  requires  that  the  density  in 
fh»  I’nirfMi  he  tore .  The  effeui.  uf  finite  density  (Holt  and  Schwartz  1963; 
Holt,  1965;  Holt,  Kawaguti  and  Sakurai),  to  first  order  is  that  of  a  pertur¬ 
bation  on  Hunter's  solution,  of  order  the  inverse  square  of  the  speed  of  the 
front,  and  is  analagous  to  the  present  work. 


2.  EQUATIONS  OF  MOTION  AMD  SIMILARITY 

The  equations  governing  the  symmetric  motion  of  a  perfect, 

inviscid,  non -heat -conducting  gas  with  constant  specific  heats  c  ,  c  , 

p  v’ 

can  be  expressed  in  characteristic  form  as 


^•(u  +  kc  )+  (u  +  c  )gj[(u  +  kc  )  =  + 


*  *  2  .* 
£  iu  c.  .  1  * 

R  7°  5R 


.  jk  .  _,Vc 

50  ,  *  50 
5T  +  u  Sr"5*  0 


where  *  denotes  a  physical  quantity. 


u  denotes  particle  velocity, 

*  *->  *  * 
c  denotes  sound  speed,  defined  by  c  «  (—^*)  *  =  > 

Sp  s  p 


s  denotes  specific  entropy 


-ix 

*V-1 


$  ,  a  measure  of  entropy,  is  defined  by  t  =  log  (- — j — ), 


Y-l  ’ 


(1) 

(2) 


and  j  =  1  for  cylindrical  symmetry 

2  for  spherical  symmetry, 
y? 

Suppose  that  U  is  the  velocity  of  a  wavefront  R  B  X(t),  moving 
into  uniform  gas.  For  the  case  of  a  strong  shock  wave  the  boundary 
values  immediately  behind  the  front,  which  are  identical  to  those  for 
a  plane  front  if  X.  is  large  in  comparison  with  the  shock  width,  are 
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2  * 
Y+l  U 


*  ^  * 

+C  rr  U 

-  Y+l 


0*  «  k  log(±  U*)  +  t* 

*•  * 

where  0^  is  the  value  of  0  at  some  reference  state  and  the  negative  sign 
* 

is  selected  if  U  is  negative.  The  assumption  that  the  shock  is  strong 

leads  to  the  neglect  of  the  undisturbed  pressure  (or  sound  speed).  Thus 

*  ★  * 

the  flow  behind  the  wave  is  determined  by  U  ,  p^,  and  since  U  has  the 

•k  \ 

dimensions  of  velocity  it  must  be  related  to  X,  t  by  U  ■  ,  where  a 

is  a  dimensionless  constant  and  t  <  0  for  the  contracting  case  (t  ■  0 
is  the  instant  at  which  the  front  is  at  R  *  0).  Hence 


so  that  the  equation  of  the  front  is 


t  m  A\ 


and  we  can  choose  A  ■  -  a  by  fixing  the  length  scale  appropriately. 


Thus  the  front  is 


*  ,1-1/a 

u  -  -  \ 


£  *  * 

Let  “  — TTn  •  ■r*ie  values  of  u  ,  c  on  the  shock,  (•  ■ 
0Ri/a 


-  1,  are 


*  _  _2_  1-1/a 

u  “  Y+l 
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then  the  equations  (5), 
d0 


di 


using  (6)  become 


expressed  in  terms  of  r,  s  after  eliminating 


2D|  §7-  (l-H-s)B,  +  (l-r-s)B 
2kD|  jJ--  (l-r+s)B+  -  (l-r-s)B_ 


(7) 


where 


D  -  (r-1) (1-r+s) (1-r-s) 

B+  ■  (r-1)  ^l-a(r+s)J  (r+ks)  +  ja(r-l)rs  +  s2 

The  equations  (7)  combine  to  give  a  single  differential  equation  for 
r  -  r(s) 

1  dr  ( 1 -r+s ) B+  +  (l-r-s)B_ 

k  ’  ds  “  (l-r-s)B  -  (l-r-s)B 
+ 


Since  the  wave  front  is  at  R  ■  0  at  the  instant  t  ■  0,  negative  values 
of  s,  which  correspond  to  negative  values  of  t,  arise  from  contracting 
fronts  and  positive  values  of  s  from  expanding  fronts. 

The  conservation  equations  across  the  front,  assumed  plane  and 
including  a  heat  release  term  are 

*  *  *  **  * 
p  (u  -U  )  «  p  (u  -  U  ) 
o  o 


where  Q 
gas  and 


*  *  *  *2  *  *,  * 
p  +  p  (u  -  U  )  -  p  +  p  (u 


o  o 


1,  *  *  2  Y 

|(a  -  U  )  +  ^ 


K-  u‘>2 


is  the  heat  release  per  unit  mass  of  gas, 
★ 

u  -  0  if  the  gas  is  initially  at  rest. 


-  uV 


+  ’  -S-  + 

Y-l  ^ 


o  denotes  the 


Q  (9) 

undisturbed 


11 


12 


linearizing  the  result  and  eliminating  derivatives  of  any  basic  terms 
by  u»  of  Luc  basic  equations .  In  terms  or  tne  dimensionless 
variables  r,  s  defined  as 

r  -  u{,  s  -  eg 


5(l-r)^  +  2(l-o)rP  -  ; 

and  hence,  in  terms  of  r(s),  s(s),  F(s), 

(1  -  r  -  s)(g.  +  k)B_  -  (1  -  r  +  .)(jjf  -  k)B+ 

(r-1)  (1  -  r  +  s)  <g  ±  k^)B±  ■  A±(gj  ±  k) 

<  r-l)||B+  -  (1-a)  (1  -  i  -  s)(^  +  k){2r(r-l)F  +  krj 


(12) 


(13) 


3.  HOMENTROPIC  SOLUTIONS 

In  the  special  case  a  ■  1,  corresponding  to  uniformly  expanding 
or  contracting  waves,  the  equation  (8)  reduces  to 


The  integral  curves  of  this  equation  are  given  in  Figure  1  (Courant  and 
Friedrichs,  1948,  page  426),  the  direction  being  of  that  of  increasing 
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time.  The  equation  has  six  singular  points 


(0,0),  (1,0),  (0.+1), 


j  I  I*  rj.  J  TRTl 


and  it  can  be  shown  that  the  nature  of  these  singularities  does  not 
depend  on  the  value  of  Y  or  whether  J  ■  1  or  2.  A  point  in  the  r,s 
plane  corresponds  to  a  path  In  the  R-t  plane.  The  possible  changes 
across  a  detonation  or  shock  front  form  a  locus  in  the  r,s  plane.  From 
the  conservation  equations  across  a  detonation  front,  with  a  constant 
heat  release  Q  per  unit  mass,  the  following  relation  between  u  ,  c 
behind  the  front  can  be  obtained 


.**  .  („*.„*)  {ijV  +  »*  .  <y-i>-3j] 

which,  in  terms  of  r,s,  becomes 


s2  -  (l-r)^*~ r+1  -  -^—^7. 

rU* 

This  is  the  equation  of  the  locus  of  the  possible  transitions  across 
a  detonation  front.  The  corresponding  shock  locus  is  obtained  by 
setting  Q  ■  0 

s2  -  (l-r)^jW^ 

which  is  an  ellipse.  In  each  of  these  equations  the  value  of  r  has  to 

2 

be  not  greater  than  which  corresponds  to  an  infinitely  strong  front. 

2  * 

The  two  curves  intersect  at  S+  where  r  ■  as  Q  is  negligible  if  U 
tends  to  infinity.  The  lines  r  ■  1  +  s  are  sonic  lines  and  are  critical 
in  that  the  direction  of  integral  curves  changes  on  crossing  them.  Thus 
no  physical  solutions  can  cross  r  ■  1  +  s.  The  detonation  locus  inter¬ 
sects  these  lines  at  D+,  which  are  the  Chapman -Jouguet  detonation  points. 
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In  the  expanding  case,  s  >  0,  an  Integral  curve  runs  from  D+  to  the 
point  (0,1),  which  corresponds  to  a  state  of  rest.  This  curve  represent* 
the  solution  given  by  Taylor's  expanding,  Chapman-Jouguet  detonation 
wave.  However,  no  integral  curve  can  be  extended  from  the  point  D  , 
corresponding  to  a  contracting  Chapman-Jouguet  detonation  front.  The 
arc  D  S  of  the  detonation  locus  corresponds  to  overdriven  fronts  and 
integral  curves  intersecting  this  arc  all  run  into  the  critical  line 
r  ■  1  +  s.  Hence  there  exist  no  uniformly  contracting  detonation  fronts, 
either  Chapman-Jouguet  or  overdriven. 

4.  THE  LIMITING  CHARACTERISTIC 

The  boundary  values  at  the  front  for  the  basic  solution  are 
given  by  (3).  However,  the  unknown  parameter  a  appears  in  the  dif¬ 
ferential  equations,  so  that  an  extra  condition  remains  to  be  found. 

This  is  obtained  by  examining  the  lines  on  which  the  solution  of 
equations  (5),  (6)  may  be  singular.  There  are  four  such  lines 

1  -  (u+c)|  -  0 
1  -  u|  -  0 


|  -  oo 


In  the  R  -  t  plane  these  lines  have  equations  of  the  form  5  *  constant 
so  that 


dR 

dt 


■ 


on  them.  Hence  the  first  pair  is 


dR 

dt 


*  * 
u  +c 


15 


i.e.,  the  positive  and  negative  characteristics  through  R  ■  0,  t  ■  0. 

Or.  i  1  >  1  ..  1  .  .  ft  IB  A  I  A  kU*  m1  n  M  It  k-lt  aU  4*  1*  A  A  M  *t  A  A  A  A  j4 

uiiuiiai  x  —  v*^  v  x x *»v-  f-Mk  ruw“  wo*  v«e..  '-**—■  vk  *6**»  -**- 

|  -  oe  is  R  ■  0,  t  $  0.  In  the  region  t  <  0  there  is  a  limiting 

negative  characteristic  (l.n.c.)  traveling  behind  the  shock  and 
reaching  R  *  0  at  the  same  instant,  t  “  0,  as  the  shock.  For  an 
arbitrary  choice  of  a  the  solution  will  be  singular  on  this  line.  Such 
a  singularity  could  exist  only  if  it  were  produced  during  the  initiation 
of  the  shock  and  precisely  on  this  limiting  characteristic.  For  this 
reason  we  shall  exclude  the  possibility  of  a  singularity  of  this  type 
and  require  the  solution  to  be  regular  on  thel.n.c.  Let  the  equation 
of  the  l.n.c.  be  £  -  ^(0  >  |,  >  -  1)  in  the  basic  flow,  and 

|  -  g1(U8R_2+2/a) 

where  6  is  a  constant,  in  the  perturbed  flow.  Thus  the  boundary  values 
^ 

of  u  ,  c  ,  0  ,  in  the  form  (11)  on  the  l.n.c.  are 


* 

u 


* 

c 


0 


* 


-  u(g1)Rl“1/a  +^)  6j8  +  uCgj)]  R”1+1/a 

^  *1 

-  c<61)»1"1/°+  6X6  +  £<6^  R'1+1/a 

k(l-i)  log  ItHKtj)  +  ftj.8  +  R'2+2/a  . 

^  li 


Also,  on  this  line 


dR  *  * 

dT"  U  ’c 


and,  from  its  equation,  we  have 


dR  _  1  C  Dl-l/a 

dT  lit* 


(3-2a)R'1+l/a  \ 
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1-1  fOL  - 14*1  fd. 

on  it.  Equating  the  coefficients  of  R  ,  R  in  these  expressions 

c  dR 

for  gives 


"  1 

5i6^ir<u"c)l  +  +  <3-2a>&  ■  ° 


where  *  ud^).  If  we  denote  by  the  value  of  s  on  the  l.n.c.  and 


r  -  r(s  ) 
o  o 


_  ,drN 

ri  "  (dS° 

S"S 


then  (14)  can  be  written  as 


r  -s  «1 
o  o 


r  -s  +  2(a-l)5  +  -4 -  5 

°  0  2s  (r.+k) 

o  1 


where  derivatives  have  been  eliminated  using  the  basic  equations. 

Since  this  line  is  a  negative  characteristic  the  variables 
there  must  satisfy  the  characteristic  condition 


d(u*-kc*> .  dt  -  £  a 


which,  on  setting  the  leading  two  coefficients  zero,  gives 


(l-i)(r-ks)  -  jr  s  - 

Ctoo  ooyCzo 


(17) 
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1 


s 

o 


6B 


erf 


2s;(ri+k) 


a-1 


(r  +s  )  +  1 
o  o 


11)21) 

y(y-D 


} 


+  8s  )  2ja-l  +  -  +  —  +  rT 

ol  J  Y  Y  1-a  o  > 

+  5  ha-3  -  -^^.7+  h  F  -o. 

L  Y  >  Y  o  o 

The  conditions  (17)  (18)  could  have  been  derived  directly  from  the 

J  „  i 

differential  equations.  We  require  — .  —  to  be  finite  on  |  *  E  . 

d|  d|  1' 

i.e.,  r  ■  1  +  $.  Hence  we  require  B  ■  0  for  a  ■  8q,  which  is  equiva¬ 
lent  to  (17).  Similarly,  the  condition  that  be  finite  for 

8  "  So  meana  that  A_  muat  be  zero  t0  order  s-so,  which  can  be  shown  to 
be  equivalent  to  (18). 


5.  THE  BOUNDARY  CONDITIONS 

•jf  tfc 

The  boundary  values  of  u  ,  c  ,  0  at  the  front  are  given  by  (10), 
from  which  we  can  deduce  the  boundary  values  of  r,  s,  r,  s,  F  there, 
taking  into  account  the  displacement  of  the  front  from  the  shock  path. 

The  equation  of  the  front  is 


5 


1  4. _ fL.  >  -2+2/a 

3 -2a  A 


Hence 


,  on  the  front 

*'‘I/0  +  ids<aF> 

0*  -  k(l-  i)  log  X.  +  0(-l)  +  ^F(-1) 


+  u(-l)l 

■-1 


5--1 


+  c(-l) 


1 


+  -*-(d*) 


X-1+1/o 


]  x-M'a 
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and  equating  these  to  (10)  and  expressing  the  result  in  terms  of  r,  8, 
r,  s,  F ,  gives 


_£)  m  _  5  2Jgy.  +  2ct-l”l  -  K 

E)  (y+l)(3-2a)  1  Y+l  i 

s(-B)  -  3~||  ^  +  QH-k-ka^  -  E1 

F(.E)  -  Ho+kp{l+S|ig>} 


The  above  boundary  values,  together  with  the  regularity  conditions 
(16),  (17),  (18)  serve  to  determine  the  solution.  The  basic  solution  for 
r  ■  r(s)  has  to  satisfy  the  differential  equation  (8),  which  contains 

2 

the  unknown  parameter  a,  and  the  boundary  values  r(sQ)  ■  r^,  r(-E)  ■  — p 

where  ro,SQ  a*e  given  in  terms  of  a  by  (15),  (17).  On  substituting  (15), 
i.e.  *o  “  1  +  s^,  into  (17)  a  quadratic  in  s^  is  obtained.  Consider  the 
expansion  for  r(s)  about  the  l.n.c. 

r(s  )  ■  r  +  r^s-s  )  +  r„(s-s  )2  +  .  .  +  r  (s-s  )n  +  .  .  . 
o  o  i  o  2o  no 


The  solution  can  be  developed  theoretically  by  substituting  this  series 
into  the  equation  (8)  and  hence  evaluating  the  coefficients  r^  by 
equating  the  coefficients  of  (s-s  )n  to  zero.  Let 


B  .  +  B  (s  -s  )  + 
0+1+0 


+  B  (s-s  )n  + 
n+  o 
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Then  the  equations  for  are 

B  »  0 

o  - 

2so(ri+  k)Bx  +  (l-r^O^-  -  0 

-2s  (r.+  k)B  -  2s  nr  B  +  r  (r,  -  k)B  -  (l-r,)nr  B 
o  1  n-  o  n  1-  n  1  o+  1  n  o+ 

+  terms  in  r  r  „  etc.  **  0,  n  >  2 

n-lJ  n-2  '  *** 

The  first  two  of  these,  which  determine  rQ,  r^,  are  quadratic  equations 

but  all  of  the  succeeding  ones  are  linear  since  B  contains  r  ,  r  , 

n-  n'  n-1 

etc.  and  is  linear  in  rfl.  Thus,  for  a  given  value  of  a,  there  are 
four  solutions  and  we  require  one  of  these  solutions  to  pass  through 
the  shock  point  for  some  particular  value  of  a. 

The  perturbation  terms  r,  s,  F  have  to  satisfy  the  differential 
equations  (13)  together  with  the  boundary  conditions  (19)  containing 
the  unknown  parameter  3,  which  measures  the  displacement  of  the  front, 
and  the  boundary  conditions  (16) (18)  containing  the  unknown  l.n.c.  dis¬ 
placement  8.  The  latter  contains  F  ,  which  is  the  boundary  value  of 
F  on  the  l.n.c.  and  is  the  third  unknown.  Thus  there  are  six  boundary  ■■n 

conditions  containing  three  unknown  parameters  8,  &,  and,  since  the 
boundary  conditions  are  linear,  the  solution  for  r,  s,  F  is  uniquely 
determined  in  terms  of  r(s). 

6.  THE  NUMERICAL  SOLUTION 

In  order  to  evaluate  the  perturbations  it  will  first  be  necessary 
to  find  the  correct  value  of  a  and  tabulate  the  basic  solution  r  *»  r(s) 
for  the  particular  choice  of  y  and  j.  We  can  avoid  the  possible  difficulty 


20 


of  —  being  indeterminate  at  a  ■  on  direct  substitution  into  the 
differential  equation  (8)  by  making  use  of  the  power  series  expansion 
for  r(s)  at  a  •  sq,  the  existence  of  which  is  ensured  by  the  regularity 
assumption.  However,  direct  computation  of  the  coefficients  is  out 
of  the  question  because  of  the  rapidly  increasing  complexity  of  the 

V 

form  of  the  equations  for  r  ,  and  each  r^  has  to  be  dealt  with 

separately.  For  this  reason  the  following  iterative  method  is  employed. 

Let  R  ,  R'  be  tabulated  functions  which  represent  r,  r'  respectively 
n  n 

as  far  as  the  term  involving  in  the  form 

R'  -  r.+2r,(s-s  )  +  .  .  +  nr  (s-s  )n_1  +  (n+l)e  (s-s  )n  +  0((a-s  )n+1) 
n  1  i  o  no  n+i  o  o 

R  “  r  +r  (s-s  >  +  .  .  +  r  <a-s  )n  +€  ..(s-8  )n+1  +  0((s-s  )n+2) 
nolo  no  n+f  o  o 


where  ^ 


is  constant. 


From  these  we  can  deduce  the  next  approximation  RW  Vi  by  substituting 
the  former  into  equation  (8),  written  as  f(r,  r',  s)  ■  0. 

Then 


2S2f 


f(R  R'  s)  -  (R  -  •)  (R* -  r')  §Jr  +  U* -r>  2 

n  n  n  Sr  n  dr  2  n  ^2 


+  <V  r><^-  r’> 


.  d2f  „  ,df  v 

where  — 2  -  0,  (^x> 

or  o 


0, 


which  we  can  write  as 


f<">  ■  Vl>  +  <s-c>”+1  + 
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In  neglecting  the  term  of  order  (s  -  s  )n+^  in  the  above  we  obtain  a 

o 

formula  for  e  ,,-r  . ,  and  hence  the  following  iterative  formula  for 

I1T  J.  irrA 


n+1 


i  -  R'  - 
n+1  n 


.  - 

/Of  \  i  .  i  \  /Of 


(s'8o>i€> 

o  1  ■» 


(20) 


The  error  coefficient  en+2  ^+2  80  0*5ta^ne<^  is  independent  of 

and  is  a  function  of  n  and  the  partial  derivatives  of  f  at  s  =  sq. 

From  (20)  we  can  tabulate  R^+1  throughout  the  range  sq  to  -E,  and  to 

obtain  R  ,,  we  use  the  following  integration  formula 
tl*T*  JL 

(■i 

j  r'd8  "  ~24  L"  rl-2  +  13rl-l  +  13rl  ‘  ri+l\  <21) 

Si-1 

where  h  is  the  step  width  and  r^  *  r'(s^).  The  formula  uses  only  points 
of  the  subdivision,  in  which  it  is  symmetric,  and  requires  one  extrapola¬ 
tion  at  each  end  of  the  range  in  the  table  of  the  derivative.  The 
relative  truncation  error  is 


7ir  where  si-2  <  aJ  <  8i+l  • 

as 

sj 

In  practice  the  total  range  from  sq  to  -E  is  roughly  0.2  so  that  only 

5  subdivisions  of  the  range  are  sufficient  to  ensure  that  the  i’-  -v 

fl“f*l 

tion  does  not  introduce  errors  of  order  (s  -  s  )  (otherwis 

o 

iteration  would  fail  to  converge  to  the  solution) .  The  extrapc  :  ■  • 

forumla  corresponding  to  (21)  is 


li+l 


4r|  -  6r^_j+  4r^_2 


-  r 


i-3 
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The  initial  approximations  are  taken  as  R|  *  r^,  Ri  ”  rc  +  rl^8"So^  an<* 
the  iteration  can  be  continued  indefinitely.  Being  iterative  the  method 
is  suited  for  programming  on  a  computer,  it  is  self -checking  to  the 
extent  that  it  can  converge  only  if  the  successive  approximations  R^, 
satisfy  the  differential  equation  more  accurately  at  each  stage.  The 
solution  is  developed  in  tabular  form,  ready  for  use  in  the  solution  for 
r,  s,  F.  Although  the  series  expansion  about  one  fixed  point,  s  *  s^, 
is  employed  the  convergence  is  very  rapid  since  the  total  range  ~  0.2. 

The  method  can  be  extended  to  the  solution  of  the  three 
simultaneous  equations  (13)  for  r,  s,  F.  Having  selected  &,  Fq 
arbitarily  we  find  r^,  sq  from  (16)and  (18).  To  form  the  initial 
approximations  to  r,  s,  F  we  require  their  derivatives  at  s  ■  s  . 

These  are  obtained  by  equating  to  zero  the  appropriate  coefficient 
in  the  expansions  of  the  differential  equations.  These  yield 


rl+  ks  ^  ■ 


-Ao+(rl+k) 
2s2B  . 

O  0+ 


ks. 


A2_(r1-k) 


B  .F.  -  -  2(l-a)<r1+  k)(2r  s  F  +  kr) 

0+  1  1  O  0  O  O 

where  A2_  is  linear  in  r^  s^  . 

Thus  we  can  tabulate  R.!  ■  r, ,  R.  *  r  +  r,  (s  -  s  )  etc.  The  iterative 

l  1  l  o  1  o 

formulae  for  R^+p  etc.,  are  obtained  by  substituting  the  n  approxi¬ 
mations  into  the  governing  equations  and  retaining  only  the  first  term, 
which  gives  three  simultaneous,  linear,  algebraic  equations  for  the  1 
corrections  -  R^,  etc.  Solving  the  equations  we  obtain  the  required 
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iterative  formulae,  and,  for  example 


r, . 


K+i  “ ft:-  pU®  +  ("+D(S.>  JL(n) 


>[{#, 


a*'  2i 


(n+1) (g~r)  M(n)(a-s0)‘2 
o 


+  - 


(sH  } 


where  equations  (13)  are  denoted  by  L,  M,  N  respectively, 


and  ,  for  example,  is  the  coefficient  of  (a  -  sq)  in  the  expansion 

of  about  s  ■  sq,  and  is  also  the  first  non-vaniahing  coefficient. 

7.  THE  BASIC  SHOCK  WAVE  SOLUTION 

For  a  given  choice  of  y,J  we  wish  to  calculate  the  appropriate 
value  of  a  and  tabulate  r(a)  from  the  l.n.c.  to  the  front.  It  remains 
to  be  settled  which  of  the  four  solutions  can  be  made  to  satisfy  the 
conditions  of  the  problem.  The  six  cases  y  ■  1.2,  1.4,  5/3  with  j  ■  1, 

2  were  computed  by  Butler  (1954).  The  same  solution  "branch"  is  taken 
in  each  of  these  cases.  In  extending  these  results  to  the  case  y  *  3, 
corresponding  to  the  motion  of  the  products  of  a  detonation,  it  is  found 
that  a  different  choice  of  branch  is  necessary.  For  this  reason  it  was 
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thought  necessary  to  examine  the  behavior  of  the  integral  curves  of  the 
dif  fcrcr.ticl  equation  with  a  view  to  examining  the  nature  of  the  change¬ 
over  and  also  the  existence  and  uniqueness  of  the  solution,  particularly 
in  the  region  of  the  changeover.  The  integral  curves  for  the  case  y  m  1.4, 
j  *  2  are  given  by  Guderley.  The  two  cases  selected  here  are  y  “  5/3, 
j  ■  2  and  y  ■  3,  J  -  2. 

The  equation  (8)  for  r(s)  can  be  written 

l  dr  r-1  _ r(r-l) (ar-1)  4-  s2^(l-a)  -  a(j+l)r^ _ 

k  ds  s  (r-i)2  -k+Cfr  (J+k+1)^  -  r(r-l)  <Or-l>  +  ks2  +  C*(l-r)} 

which  has  nine  singular  points.  Thure  are  three  on  the  r-axiB  P^(0,0), 
P^(0,1)  and  (0,^)  and  three  in  the  region  s  <  0 


W-  ' +  V> 


r A,  ,  i  +  .  ) 

J  O-  O  - 


VS'  a(j+k+l)  * 


using  Guderley' s  suffices.  The  remaining  three  are  the  mirror  images  of 

P2'  P3*  P5  in  t*ie  r"ax:''8  and  correspond  to  expanding  flows.  The  quantities 

s  ,  are  the  two  roots  of  the  quadratic  for  s  and  S  is  the  negative  solu- 
0+  o 

tion  of 

S2  -  r(r-l)(ar-l) 
k  + 

where  r  "  a(jfe-TT  ’ 

The  behavior  of  the  integral  curves  is  found  by  determining  the 
nature  of  these  singularities,  the  region  of  interest  being  s  <  0, 


O 
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0  <  r  <  1.  In  these  calculations  the  correct  value  of  a  was  used.  The 

~ _ —  -  —  _  —  1 _ i _ L  .  J  J  _  T5  J _  O  f _  >1  .  .  _  ..  __  m  i  r\  ,  _  n  1  «  « 

VCO  Bfc  V  or.ct.vm;u  Ait  a  i^ulc  a  IUI  UIC  V.  a  O  C  J  —  _/  /  J  ^  j  —  A.  nu  curves 

change  direction  on  crossing  the  line  r  ■  1  +  s,  except  for  the  two 

limiting  ones  through  each  of  P2  and  P^,  which  represent  the  four  solu¬ 
tions  which  are  regular  on  the  l.n.c.  We  require  a  curve  which  starts 
at  the  shock  point  and  passes  through  P2  or  P^  and  also  through  the 
origin  P^,  which  corresponds  to  t  ■  0.  On  this  curve  time  must  increase 
from  the  shock  point  to  P^.  From  the  sketch  it  is  seen  that  there  are 
two  such  curves,  one  through  P2  and  the  other  through  P.^.  One  of  these 
has  to  be  made  to  pasB  through  the  shock  point  for  some  choice  of  at. 

For  values  of  y  in  this  neighborhood  it  was  found  in  practice  that  an 
appropriate  solution  was  found  by  selecting  the  curve  through  P^  and 
the  curve  through  P^  could  not  be  made  to  pass  through  the  shock  point. 

The  sketch  of  the  integral  curves  in  the  case  y  *  3,  j  ■  2  are 
given  in  Figure  3.  In  performing  the  calculations  to  determine  the 
nature  of  the  singularities  the  correct  value  of  a  was  used.  In  thiB 
case  P2  and  P^  are  both  nodes  and  P^  is  a  saddle  point,  below  the  line 
r  ■  1  +  s.  The  two  curves  running  towards  P,.  can  be  discounted  but 
there  is  no  obvious  choice  between  the  remaining  two.  Again  only  one 
solution  was  found,  the  curve  through  P^  being  selected  in  this  case. 

Investigation  of  the  Integral  curves  does  not  settle  the 
issue  of  uniqueness  of  the  solution.  For  given  y,j  there  is  a  range  of 
values  <  a  <  a2  for  which  sq  is  imaginary.  The  range  0  <  a  < 

never  yields  any  solutions.  For  j  ■  2  as  y  is  increased  from  1.2  the 
correct  value  of  a  approaches  a 2  and  the  correct  values  of  sq  approach 
each  other.  For  some  critical  value  of  y,  y^  say,  the  roots  are  equal 
and  the  transition  from  one  branch  to  the  other  occurs  at  y£.  For  any 
given  value  of  y,  a2  is  that  value  of  a  for  which  P2,  P^  coincide.  For 
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values  of  a  <  these  two  singular  points  are  complex  and  so  no  regular 

solutions  can  oe  continued  across  the  line  r  ■  1  +  s  to  the  origin.  As 

a  is  increased  from  P^  and  P^  separate  and  move  along  the  segment  of 

r  ■  1  +  8  as  far  as  (0,1),  (-1,  0)  when  a  ■  1. 

We  require  a  solution  through  either  P^  or  P^,  the  solution  and 

the  positions  of  the  two  points  depending  on  the  value  of  a,  and  also 

through  the  shock  point,  the  position  of  which  depends  on  y  only.  In 

Table  1  d(P2)  denotes  the  discrepancy,  for  the  spherical  case,  between 

the  solution  obtained  by  Integrating  from  P2  as  far  as  s  *  -E  and  the 
2 

required  value  of  there.  For  Y  ■  1.865  d(P3)  has  a  zero  in  the  given 
range  and  this  zero  corresponds  to  the  actual  solution.  Apparently  d(Pg) 
has  no  zero.  The  situation  is  reversed  in  the  case  Y  “  1.875,  the  point 
P2  being  appropriate  in  this  case.  Thus,  1.875  >  Yc  >  1.865.  The  tran¬ 
sition  at  y  *  Vc  takes  place  smoothly  and  there  i ,p  no  apparent  physical 
significance  to  the  case  Y  ■  Yc* 

For  a  given  y  the  rootp  for  sq  are  monotonic  in  a  in  the  range 
0  <  a  <  1,  so  that  P2,  P3  vary  continuously,  without  repetition,  along 
the  arc  r  ■  1  +  s  as  a  varies  between  a2  and  1.  Together  with  the  results 
of  Table  1  this  suggests  the  following  behavior.  For  a  m  a^,  P2  and  P^ 
coincide  and  the  single  integral  curve  through  them  separates  the  area 
0  <  r  <  1,  r  >  1  +  s  into  two  distinct  regions.  As  a  is  increased  the 
two  integral  curves  through  P2,  P^  must  lie  wholly  within  each  of  these 
regions  so  that  points  in  the  lower  region  may  be  reached  from  P^  for  some 
value  of  a  and  those  above  from  P2 .  It  seems  likely  that  no  two  curves 
through  one  of  P2,  F.  will  intersect  for  distinct  choices  of  a.  If  this  is 
so  then  the  solution  will  be  unique  for  all  y.  The  choice  between  P2,  P^  is 
determined  by  whether  the  shock  point  lies  above  the  limiting  integral  curve 
through  the  point  formed  by  the  merging  of  P2,  Pj.  Apparently  for  Y  <  Yc 
the  shock  point  lies  below  this  curve,  and  above  it  for  Y  >  Yc 
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The  results  for  the  eight  cases  V  ■  1.2,  1.4,  5/3,  3  with  j  ■  1,2 
are  given  in  Table  2  along  with  those  given  by  Whitham's  approximate 
method,  to  be  described  later,  for  comparison.  The  case  y  m  3  has  not 
been  studied  previously. 

8.  THE  PERTURBATION  SOLUTION 

The  solution  for  the  perturbation  s  are  now  obtained  by  integrating 
the  three  simultaneous  equations  (13)  for  r,  s,  F,  subject  to  the  boundary 
conditions  at  the  front  and  the  l.n.c.  The  function  r(a)  and  the  para¬ 
meter  a  appearing  in  (13)  are  now  known.  As  for  the  basic  solution  we 
develop  the  solution  away  from  the  l.n.c.,  having  satisfied  the  regularity 
condition  there,  as  far  as  the  front.  To  do  so  wc  select  arbitrary  values 
of  8,  F  ,  which  determines  r,  s,  F  at  s  ■  a  ,  and  continue  the  solution 

O  O' 

to  b  ■  -E,  where  the  conditions  will,  in  general,  not  be  satisfied  by 
the  present  solution.  Suppose  we  have  found  two  such  linearly  independent 
solutions,  corresponding  to  choices  6^°^,  F^°^  and  8^,  for  the 

values  of  8,  F  .  let  denote  the  values  of  r  at  s  *  s  , -E 

respectively  of  the  o  solution.  The  boundary  conditions  at  the  front, 
given  by  (19)  can  be  written  as 

rH  -  Al  P  -  K 

SH  "  A2  P  'E' 

F„  *  A  p  +  H  ,  where  P  is  unknown. 

H  J  O 

Let  us  take  a  linear  combination  of  the  two  numerical  solutions  and 
satisfy  the  above  conditions.  Thus 
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x^o)  +  y;^1>  -  Ajp  -  k 

xI<o)  +  ys^1)  »  a23  -  E'  (22) 

XF„°^  +  YF^  ■  A_3  +  H 
n  n  J  0 

which  can  be  readily  solved  by  X,  Y,  |3  so  that  we  can  thus  evaluate  p. 

The  solution  for  r,  s,  F  appropriate  to  the  boundary  conditions  can  be 

obtained  by  performing  the  integration  from  the  l.n.c.  with  the  correct 

values  of  8.  F  which  are  given  by 
o 

F  -  XF^  +  YF^ 

0  0  o 

X(r(o)+  s^)  +  Y(r(1^+  3(1)) 

R  -  o  o  o  o 

2(a-l)  +  (rl~1)Bo+ 

2  2Q(v1+k) 

\ 

The  quantities  K,  E',  H  appearing  in  (22)  depend  upon  the 
2  0 

values  of  Q,  cq  .  However,  the  dependence  1b  linear  so  that  all  that 

is  necessary  is  to  evaluate  two  solutions  due  to  linearly  independent 

choices  of  Q,  cq  ,  For  simplicity  we  can  take  Q  *  .1*  c  *  0  and 

Q  “  0,  cq  *1,  the  former  corresponding  to  a  detonation  front  and 

the  latter  to  the  correction  due  to  counter-pressure.  The  solution  in 

a  specific  case,  due  to  either  or  both  of  these  effects,  is  found  by 

taking  the  appropriate  combination  of  these  two  solutions.  For  y  ■  1.2, 

1.4,  5/3,  3  and  j  ■  1,2  we  have  sixteen  distinct  cases.  The  results 

for  these  are  given  in  Tables  3,  4,  5,  and  6  along  with  those  obtained 

ft  ft 

by  the  approximate  method  of  Whitham.  The  boundary  values  of  u  ,  c  at 
the  front  are  given  by 
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u*  -  -  ^  X1-1/“{1  +  <(.-  *!•  K>d+l/“] 
c"  -  EX*-1,C,|  1  +  (*■  §V2+2/“] 

The  coefficients  f3  -  K,  (3  +  |—  are  also  tabulated. 

9.  THE  RESULTS  AND  THE  WHITHAM  SIMPLIFIED  ANALYSIS 

Before  discussing  the  results  It  will  be  of  interest  to 
evaluate  the  solution  by  the  approximate  method  in  the  form  given  by 
Whitham.  It  is  known  that  this  gives  remarkable  accuracy  in  estimating 
a  (Chisnell,  1957 j  Whitham,  1958).  Chisnell  employed  his  "shock-area" 
rule,  which  he  formulated  for  shod;  waves  in  channels  of  slowly  varying 
cross-section,  in  the  evaluation  of  a  for  y  ■  1.2,  1.4,  5/3  and  j  «  1,2 
and  compared  the  results  with  Butler's.  Whitham  obtained  Chisnell' s 
results  by  assuming  that  the  characteristic  conditions  to  be  satisfied 
behind  the  shock  will  be  satisfied  by  the  boundary  values  there.  This 
method  will  be  applied  to  the  present  problem. 

The  characteristic  condition  to  be  satisfied  behind  the  front 
is 

*  *  *  *  -1  i  *  .* 

d(u  -  kc  )  «  ju  c  R  dt  -  dt) 


and  the  boundary  values  there  are 

* 


r~r  R1"1^01  +  (K- 
Y+l 


Mr)  R'1+1/a 

Y+l 


* 

c 


ER1”1^  +  <E'+  pE>  R’1+1/a 


0*  -  k(l-  £)  log  R  +  (H  +  kp)  R"2+2/a 
u  o 
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On  substituting  these  values  into  the  characteristic  condition, 


_  J  _  ^  «.  t _ C _ 1.1. 

auu  uull  ije,  LUC  1QLL  Lilac 


dR 

dt 


* 


* 

u  ,  we  ifuLdxii  a  puiy iicnua.ax 


„-2+2/a 

IV  , 


Setting  the  first  term  zero  gives  Whitham' 8  formula  for  at 


1 

a 


1 


_JE 

Y+l _ 


(23) 


Equating  the  coefficient  of  the  second  term  to  zero  yields 

h±nL  -  4. 


e2k  - 


^EH  +  K  -  —  E' 


P  - 


-A-E(-i—  +  E) 

Y+l  Y+l 


JOL 


4(Y+1  +  YJ 


(24) 


These  approximate  results  neglect  the  effect  of  disturbances  reaching  the 

shock  from  behind,  due  to  the  characteristic  condition  not  being  applied 
#• 

correctly.  The  changing  surface  area  of  the  shock  is  accounted  for.  The 
area  of  the  front  is  proportional  to  R^  which  results  in  the  exponent  in 
the  power  law  for  the  shock  speed,  i.e.,  1-^,  being  proportional  to  j. 

The  perturbation  solution  for  p,  which  arises  from  energy  terms  pro¬ 
portional  to  volume  and  Independent  of  the  geometry  of  the  syBtem,  is 
independent  of  j.  The  result  (23)  for  a  is  very  accurate  because  the 
propagation  of  the  shock  is  largely  governed  by  the  focussing  effect, 
due  to  its  surface  area  diminishing,  and  is  affected  little  by  other 
disturbances.  This  is  not  the  case  for  the  perturbations  neither  of  which 
(heat  release  and  initial  pressure)  are  geometric  effects  and  the  results 
for  0  are  much  less  accurate  than  those  for  at.  A  graph  of  the  approxi¬ 
mate  results  for  0  is  given  in  Figures  4  and  5. 

From  the  results  obtained  by  the  full  analysis  it  is  seen  that 
for  given  j,  Q,  cq  there  is  always  a  change  in  sign  in  p,  considered 
as  a  function  of  y.  Thus  the  introduction  of  either  of  the  two  effects 
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can  produce  an  Increase  or  a  decrease  in  the  speed  of  the  front, 

depending  upun  Luc  value  of  y-  In  each  case  p  ■  o  tor  some  value  of 

Y  and  in  this  case  the  perturbation  of  the  front  speed  in  of  order 

X  The  critical  values  of  y  were  found  by  the  approximate  methc^ 

2  2 

to  be  1.30  for  Q  ■  0,  cq  b1  and  1.24  for  Q  ■  1,  cq  ■  0.  The  graph 
of  0  from  the  approximate  method  follows  the  behavior  of  the  correct 
values  fairly  closely,  and  the  former  would  appear  to  be  sufficiently 
accurate  to  estimate  the  critical  values  of  y  to  two  decimal  places. 

If  the  initially  uniform  medium  considered  so  far  is  replaced 

k _ _  ju 

by  a  medium  initially  at  rest  but  having  variable  density,  pq  CC  R  , 

*  -m/2 

say,  where  m  is  a  constant,  then  the  initial  Bound  speed  c^CCR  . 

* 

The  similarity  hypothesis  will  still  hold  provided  cq  remains  small 

.tff  m  i 

relative  to  U  ,  i,e.  j  <  —  -  1,  and  the  solution  is 

U*  -  -  x1_1/q!ci  +  pX-1+1/'a-m/2) 


where  the  value  of  0  is  identical  to  that  in  the  case  cq  «  1,  Q  ■  0, 
computed  previously.  The  coefficients  of  the  perturbations  do  not 
differ  from  the  previous  solution,  the  only  difference  being  in  the 
power  of  X. 

The  author  gratefully  acknowledges  the  guidance  of  Mr.  D.  S. 
Butler  throughout  the  course  of  this  work  and  in  particular  for 
suggesting  the  numerical  solution.  The  computations  were  performed 


on  the  Atlas  Computer  of  the  Science  Research  Council,  Didcot, 


Bekshire . 
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Table  1 


Y  -  1.865 

Y  -  1.875 

a 

d(P2) 

d(P3) 

a 

d(P2) 

d(P3) 

0.674453 

so  lmag. 

so  lmag. 

0.6738558 

so  lmag. 

ao  lmag. 

0.6744535 

-  0.00033 

0.6738559 

+  0.00206 

+  0.0033 

0.674454 

-  0.00557 

+  0.00066 

0.6738560 

+  0.00132 

+  0.0039 

0.674456 

-  0.00798 

+  0.00308 

0.6738562 

0.674460 

-  0.0109 

0.67386 

|  -  0.00387 

1  +  0.0092 

0.675 

-  0.0749 

+  0.0751 

L 

0.675  i 

1 

-  0.1003 

+  0.118 

Table  2 


1.2 

0.861163 

0.859762 

0.757142 

0.754021 

1.4 

0.835323 

0.835373 

0.717174 

0.717288 

5/3 

0.815625 

0.816043 

0.688377 

0.688654 

3 


0.775667 


0.772661 


0.636411 


0.629542 
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selected  *o  that  no  security  classification  ia  required.  Idtntl 
flare,  auoh  ae  equipment  model  designation,  trade  name,  mllltaiy 
project  cede  name,  geographic  location,  may  be  used  aa  key 
words  but  will  ba  followed  by  on  indication  of  technical  con¬ 
test.  The  assignment  of  Ueks,  rales,  end  weight*  ia  optional. 
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